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We present a detailed theoretical study of the influence of linearly polarised intense ter- 
ahertz (THz) laser radiation on energy states of hydrogen-like impurities in semiconductors. 
^3 ;The dependence of the binding energy for ground (Is) and first excited (2s) states, E\ s and 
^ '.E 2s , on intensity and frequency of the THz radiation has been examined for a GaAs-based 
O .system. It is found that E\ s , E2 S and E 2s — E\ s decrease with increasing radiation intensity 
& or with decreasing radiation frequency, which implies that an intense THz field can enhance 
ionisation of dopants in semiconductors. Our analytical and numerical results show that one 
- of the most important results obtained by A.L.A. Fonseca et al [phys. stat. sol. (b) 186, 
g ;K57 (1994)] is incorrect. 

X 
S3 



1. Introduction 

The ionisation of dopants (i.e., donor and acceptor impurities) provides the major source 
of carriers (i.e., electrons and holes) in semiconductors. The ionised dopants are also the 
major sources of electron-impurity scattering in semiconductor devices, which determine the 
transport and optical properties of the device systems at low-temperatures. Hence, the 
investigation of ionisation of dopants (such as binding energies of donors and acceptors, 
transition energies and probabilities among different impurity states, etc.) is fundamental in 
understanding almost all physically measurable properties in semiconductors. In the absence 
of an intense electro-magnetic (EM) radiation field, the binding energies of dopants and the 
transition energies among different impurity states in popularly used semiconductor materials 
are known [1] and the theoretical approaches to calculate these impurity states are well- 
documented [2]. 

It should be noted that in semiconductor materials such as GaAs, Ge and Si, the binding 
energies of donor and acceptor impurities are of the order of terahertz (10 12 Hz or THz) pho- 
ton energies [3] so that an intense THz radiation can affect strongly the impurity states. With 
development and application of coherent, high-power, long-wavelength, frequency-tunable 
and linearly polarised radiation sources such as THz or far-infrared (FIR) free-electron lasers 
(FELs) [4], it has now become possible to measure the effect of an intense laser radiation on 
ionisation and perturbation of dopants (especially shallow impurities) in different semicon- 
ductor systems. In recent years, using THz FELs [such as Free Electron Laser for Infrared 
experiments (FELIX) in The Netherlands and CW FELs at UCSB] as intense radiation 
sources, THz-photon-induced impact ionisation in InAs heterostructures [5], time-resolved 
shallow donor spectrum in Si-doped GaAs [6] and Lyman transitions in Be-doped GaAs [7] 
have been investigated through, e.g., transport and/or photoconduction measurements. The 
results obtained experimentally indicate that in the presence of intense laser radiation such 
as FEL fields, i) impurity states in different semiconductor systems are perturbated by the 
intensity and frequency of the THz laser fields [5-7]; ii) photoconduction experiments are 
more sensitive than optical measurements for detecting transition energies of impurity states 
in semiconductors [6,7]; and iii) some interesting intense radiation phenomena, such as im- 
pact ionisation of dopants [5] and splitting and broadening of the impurity spectrum [7], can 
be observed. In order to understand these fundamentally new experimental findings and to 
predicate new radiation phenomena, it is essential to know theoretically how an intense laser 
field affects the binding energies of impurities in semiconductors and is the prime motivation 
of the present study. 

In this paper, we study binding energy of hydrogen-like impurities in bulk semiconductors 
in the presence of intense THz laser fields. For semiconductor materials such as GaAs, there is 
an absolute conduction band minimum at K = with T 6 symmetry and the electron-effective- 



mass and dielectric constant can be considered to be isotropic. Therefore, the impurity 
states in materials such as GaAs are strictly scaled-down versions of those in the hydrogen 
atom. The small electron-effective-mass coupled with the large static dielectric constant 
result in small binding energy for donors in these materials [3] . It has been noticed by atomic 
physicists in their early work that intense UV radiation can affect strongly the structure of 
atomic hydrogen [8] . Using theoretical approaches proposed and developed in Ref. [8] , Nunes 
and co-workers have investigated the influence of intense laser radiation on hydrogen-like 
impurity states in bulk semiconductors [9] and in semiconductor quantum well structures 
[10]. However, in Ref. [9], only the ground-state energy at high-intensity radiation was 
presented analytically. We find that it is hard to understand the result given by Eq. (13) in 
Ref. [9], because it shows that at high-intensity radiation limit (ie, «o —> oo) the ground-state 
binding energy Eq — > — R* (binding energy in the absence of the radiation) with R* being 
the effective Rydberg constant. In order to see how impurity states in a semiconductor are 
perturbated within the frequency and intensity range of the current generation of the THz 
FELs, we think it is necessary to re-examine the dependence of the impurity binding energy 
in a semiconductor on frequency and intensity of the THz fields. Moreover, since most of the 
experimental work in this area has been carried out using photoconduction experiments [6,7] 
which measure the transition energies among different impurity states, it is also necessary to 
know how an intense THz radiation affects the excited impurity states in a semiconductor. In 
Section 2, we will first briefly introduce the theoretical approaches used for the calculations 
then present analytical results of the binding energies for ground and first excited impurity 
states. The numerical results will be presented and discussed in Section 3 and the conclusions 
drew from this study will be summarised in Section 4. 

2. Analytical results 

In this work, we consider the situation where a laser field with a vector potential A(t) 
is applied along the x-direction of a semiconductor and the radiation field is polarised lin- 
early along the z-direction. Under the effective mass approximation, the time-dependent 
Hamiltonian to describe the electron-impurity system is given as 

Hit) = ^[Pl+P 2 y + (Pz ~ eA{t)f] - (1) 

Here, m* is the electron-effective-mass, p x = —ihd jdx is the momentum-operator along 
the x-direction, under the dipole approximation A(t) = (F /o;)sin(a;t) is the vector potential 
induced by the radiation field with uj being the radiation frequency and F the electric field 
strength of the radiation field, R = (r,z) = (x,y,z), and e 2 / kR is the Couloumb potential 
induced by electron- impurity interaction with k being the dielectric constant. It can be seen 
that due to time- dependent nature of the radiation field, in principle, we have to solve time- 
dependent Schrodinger equation: [ihd/dt — H(t)]^f(H,t) = to obtain the electronically 



bounded impurity states. However, at present, there is no simple analytical solution to this 
problem [11]. Here we employ a tractable theoretical approach proposed and developed by 
atomic physicists [8] to study the binding energy of the impurity states. The key issue of this 
approach is to find laser-dressed potential energy. It is found that if the coordinate along 
the z-axis is shifted by z' = z + zo(t) with zo(t) = — (eF /m*tv 2 )cos(tvt), the Schrodinger 
equation becomes 



dt 

where R/ = (r, z') and 



Pl+ ^ p2z ' +V(R\t)]*(R',t), (2a) 
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V(K',t) = = + 1 - WJ . (26) 

V ' K^r 2 + [z> - z (t)} 2 2m* v ' 

After averaging V(R',t) over a period of the radiation field, we obtain laser-dressed potential 
as ^ 

V(R') = -£- [ dt y(R / ,t) = J E em -^[|R / + R |- 1 + |R / -R r 1 ]. (3) 

Here E em = (eF ) 2 /4m*tv 2 is the energy shift induced by the radiation field due to dy- 
namical Franz-Keldysh effect [12] and Ro = (0, 0, eF /m*uj 2 ). It should be noted that in 
atomic physics, E em is also called ponderomotive energy [13]. If we replace V(Rf,t) in 
Eq. (2a) by T^(R') given by Eq. (3), the time-dependent problem becomes approximately a 
time-independent one and the solution is ^f(R.',t) = e - t ( E N+E em )t/n^f ^(R/) ; where the wave- 
function ^jv(R') and the energy spectrum En are determined by (H — En)^n(R') = 
with 

Ho = Pi tA±A _ |!(|R- + Bo|-' + |R' - R |-'). (4) 

In contrast to Refs. [8-10], in the derivations shown above we have included the terms 
associated with Fq which, as can be seen, only contributes to the energy shift induced by 
the radiation field alone. The main merit of this approach is that now one can convert a 
time-dependent problem into a time-independent one and, thus, one can use, e.g., variational 
method to calculate the energy states and the corresponding wavefunctions. It should be 
noted that the present theoretical approach is a generalization of those proposed by the 
atomic physicists. As has been pointed out by Ref. [13], this approach is the lowest order 
approximation to the atomic problem. Furthermore, the results given by Eqs. (3) and 
(4) have been documented by Refs. [9,10] for the case of the hydrogen-like impurities in a 
semiconductor. 

In the present work, we employ the wavefunctions for a hydrogen atom as the trial 
functions for ground state (ie, ls-state) and first excited state (ie, 2s-state) of the donor 
impurities, which reads 

*i fl (R ; ) = ae-^K (5a) 



and 

tt 2a (R/) = c(l + bR'/al)e-^ R '' 2 <\ (56) 

where (3 and 7 are variational parameters, clq = h 2 K/m*e 2 is the effective Bohr radius, 
a = (/3 3 /7rao 3 ) 1/2 is determined by < ls\ls >= 1, and 6 = -{(3 + t/2)/3 and c = 
7 5/2 /[87raQ 3 (126 2 +667+7 2 )] 1/2 are determined respectively by < ls\2s >= and < 2s\2s >= 
1. After introducing these trial functions to the variational method, i.e., En =< N\Hq\N >, 
the binding energies of the ls-state and 2s-state are obtained as 



Eis — —R y 



% - (3 2 - ^{l + al(3)e- 2 <?\ (6a) 
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R* y (4/a*)(246 2 + 126 7 + 2 7 2 - Ae~ a ^) - i 2 {4b 2 + 26 7 + 7 2 ) (au . 
E2s ~- — 1262 + 667 + 7 2 ' { ) 

where «q = eF /m*uj 2 aQ, R* = m*e 4: /2h 2 K 2 is the effective Rydberg constant, and 

A = 24b 2 + 1267 + 2 7 2 + ao7(186 2 + 867 + 7 2 ) + 2af 6 7 2 (36 + 7) + a^ 3 6 2 7 3 . 

Furthermore, the variational parameters (3 and 7 are determined respectively by 

P-0- + 2a* p)e~ 2 < f3 = and dE 2s /d~f = 0. 

The results shown above indicate that the effect of the radiation field on binding energy of 
the impurities is mainly achieved through a dimensionless factor «q = e 3 Fo/h 2 KU! 2 ~ Fq/lo 2 . 
For high-frequency and/or low- intensity radiation so that < 1, we have 

E la ~-R;(l- -a?) and E 2s ~ —jL (l - -a* 2 ) , (7) 

which suggests that the binding energy decreases quickly with increasing radiation intensity 
or with decreasing radiation frequency when ckq ~ 1. When = 0, E\ s = —Ry and 
E 2s = —R y /^ are well-known results obtained in the absence of the radiation field. For 
low- frequency and/or high-intensity radiation, entailing ckq » 1, we find that there is no 
simple analytical expression for both E\ 8 and E 2s , in contrast to Eq. (13) obtained by Ref. 
[9]. When ckq — > 00, we find (3 — > and 7 — > so that £i s — >• and E 2s — > (note that a 
zero binding energy implies that all impurities at this state are ionised). These theoretical 
results indicate that E\ s and E 2s are altered from respectively — R* and — R*/4 at zero field 
to zero at high- field limit. One of the most important results obtained by Ref. [9] is that at 
high- intensity radiation fields, the ground-state binding energy is given by Eq. (13) which 
shows E\ s — ► —R y when ao — > 00 and E\ s < —R y - Our results here suggest that Eq. (13) 
obtained by Ref. [9] is incorrect. 



3. Numerical results 

Below we present numerical results for semiconductor materials such as GaAs. The 
material parameters for GaAs taken within the calculations are the effective-electron-mass 
ratio m* /m e = 0.0665 with m e being the electron rest mass and the static dielectric constant 
k = 12.9. The dependence of binding energies E\ s and E 2s as well as transition energy 
E 2s — Eis on THz laser radiation fields is shown in Figures 1-4. From these results, we see 
that: 

a) in the presence of the radiation fields, E ls and E 2s are altered from respectively — R* 
and —R*/4: at low-field limit to zero at high- field limit and E 2s — E\ s is altered from 
3.R*/4 at low-field limit to zero at high- field limit. This confirms that Eq. (13) obtained 
by Ref. [9] is incorrect; 

b) with increasing radiation intensity and/or decreasing radiation frequency, E\ S) E 2s and 
E 2 s — Eis decrease; 

c) E 2s depends a bit weakly on the radiation field than Ei s does, which can be understood 
byEq. (7); 

d) the strong effect of the radiation field on binding energy and transition energy can be 
observed at ~ 1; and 

e) Eis, E 2s and E 2s — Ei s depend more strongly on radiation frequency than on radiation 
intensity because ~ Fq/uj 2 . 

The results discussed above indicate that in the presence of the intense THz laser fields, 
the binding energies of the impurity states can be reduced and the impurity spectrum in 
semiconductors can be shifted significantly by the radiation. 

It should be noted that the current generation of the FELs can provide intense THz 
radiation sources in the frequency and intensity range / ~ 0.1 — 10 THz and Fq ~ 0.1 — 100 
kV/cm [14] so that the condition «q = e 3 Fq / ft 2 klv 2 ~ 1 can be satisfied by most of the 
popularly used semiconductor materials. We therefore believe that it has now become possible 
to investigate the effects of the intense laser radiation on ionisation and perturbation of 
dopants in semiconductor systems by using current generation of THz or FIR FELs. 

4. Further remarks 

In the present study, our calculations have been performed for shallow-impurities in 
semiconductor materials such as GaAs, because these impurities are hydrogen- like. It is 
well known that in semiconductors, deep-impurities may not be hydrogen-like because their 
states depend strongly on lattice structure of the material systems, such as the symmetry 
of the crystal potential, electronic band structure, etc. To examine the effect of intense 
laser radiation on deep-impurity states in a semiconductor, the above mentioned theoretical 
approach has to be modified. 



To our knowledge, except those theoretical results reported by Refs. [9] and [10] for semi- 
conductor systems, so far very few theoretical and experimental results have been reported 
regarding the direct effects of the intense laser radiation on binding energies of electrons in 
atomic systems (two recent reviews are given in Refs. [13] and [15]) or of hydrogen-like im- 
purities in semiconductor systems. The main reason behind this is that when an electronic 
system (atoms or semiconductors) is subjected to intense laser fields, the investigation of the 
electronic states (or electronically bounded impurity states) is essentially a time- dependent 
problem. As a consequence, these states are no longer the eigenstates and transitions among 
these time-dependent states can occur. Therefore, it is very hard to measure experimentally 
the binding energy for a certain atomic or impurity state when an intense laser radiation is 
present. Normally, the results obtained from atomic physics experiments [13,15] are contri- 
butions from all possible states and transitions among these states. 

The analytical and numerical results shown in this paper are what one can get from 
the simplest theoretical model. As has been pointed out by Ref. [13], the approach used 
in the present study is the lowest order approximation to deal with the problem. Under 
such an approximation, the systems can be described in terms of a set of quasi-eigenstates of 
atom (or impurity) and laser field and there is no transition among these quasi-eigenstates. 
Namely, under this approximation, the system is "stable to ionisation". Thus, the results 
obtained from this model can be useful in understanding atomic physics phenomena such as 
free-free transitions and dichotomy of the atomic states [8] . Furthermore, it can be shown [13] 
that the conditions under which the approach can be applied are Too 3> E em and RqUJ S> 1. 
Therefore, the very low-frequency and high-intensity results shown in Figs. 1-4 in this paper 
may not be the case. To examine the effect of these extreme radiation conditions on states of 
hydrogen-like impurities in a semiconductor, the contributions from higher order terms (see 
Eq. (28) in Ref. [13]) have to be included. The inclusion of these higher order contributions 
requires considerably further analytical and numerical calculations and, therefore, we do not 
attempt it in the present study. 

5. Conclusions 

In this paper, we have re-examined how a linearly polarised intense laser field affects the 
binding and transition energies of shallow-donor impurities in a bulk semiconductor. The 
present study has been conducted on the basis of a theoretical approach proposed by Ref. 
[8] and documented by Refs. [9,10]. We have obtained results of binding energies for both 
ground and first excited impurity states and thus we can look into the effects of the radiation 
field on transition energy between Is and 2s states, in conjunction with experimental work 
carried out by photoconductivity measurements [6,7]. The numerical and analytical results 
obtained from this work show that Eq. (13) in Ref. [9] is incorrect. 



We have found that the strong influence of the radiation field on binding and transition 
energies of the shallow impurity states in bulk semiconductors can be observed when the 
condition «q = e 3 F /h 2 KLV 2 ~ 1 is satisfied. For most popularly used semiconductor mate- 
rials such as GaAs, Si and Ge, these radiation conditions have been realised by the current 
generation of the THz or FIR FELs. We therefore hope that the phenomena studied and 
predicated in this paper can be verified experimentally. 
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Figure captions 



Fig. 1: Binding energies, E\ 8 and -E^s, as a function of THz radiation frequency (/ = lo/2tt) 
for different radiation intensities. R* is the effective Rydberg constant and for GaAs R* = 
5.44 meV. 

Fig. 2: Transition energy between 2s and Is impurity states, i?2s — E\ s , as a function of 
radiation frequency for different radiation intensities. 

Fig. 3: Eis and E 2s as a function of radiation intensity for different radiation frequencies. 
Fig. 4: E2.S — Eis as a function of radiation intensity for different radiation frequencies. 
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